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Abstract 

The two state molecular orbital (2-MO) model of the phenyl based semiconductors is used 
to calculate the low-lying spectra of the and B^^ states of poly(para-phenylene) (PPP). The 
model parameters are determined by fitting its predictions to exact Pariser-Parr-Pople model 
calculations of benzene and biphenyl, and it is solved using the density matrix renormalisation 
group method. It is shown that there exists a band of ^B^^ ('s'-wave) excitons below the band 
states. In the long chain limit the lowest exciton is situated 3.3 eV above the ground state, 
consistent with experimental data. The calculated particle-hole separation of these excitons 
indicates that they are tightly bound, extending over only a few repeat units. The lowest band 
state is found to be a covalent 2^j4^ state, whose energy almost coincides with the charge gap 
Eg- Lying just above the 2^A'^ state is a band ^B^^ state (the n^B^^ state). The particle-hole 
separation of the band states scales linearly with oligomer size. The binding energy of the f ^-B^ 
exciton is determined rigorously as 0.74 eV. 

The dipole matrix elements and oscillator strengths for the transitions between the lowest 
^A^ and ^-B^ states are calculated and the NLO properties of PPP, such as electroabsorp- 
tion (EA) and third harmonic generation, are investigated. A comparison of the EA spectrum 
with the experimental data show that the main features of the experimental spectrum are well 
described in the 2-MO Hamiltonian. 

Only five states account for most of the calculated EA. These are the 1^A+, l^Bi" , 2M+, 
n^B^^ and another band ^A^ state, the k^A'^, thus confirming the essential states model. An 
analysis of the particle excitation weight of these states indicates that they are predominately 
single particle in character. 

PACS numbers: 42.70Jk, 71.20Rv, 78.66Qn, 71.35Cc 



1 Introduction 

Since the first light-emitting device based on poly(para-phenylenevinylene) (PPV) was reported 
the non-hnear optical (NLO) properties of conjugated polymers have been extensively investigated. 
Amongst the numerous systems studied, poly(para-phenylene) (PPP), being a linear chain of phenyl 
rings, possesses one of the simplest structures. However, its electronic structure and the nature of 
the blue light emission arc still controversial. First-principles local-density approximation studies 
by Ambrosch-Draxl et al. |^ suggest that the optical properties of PPP can be explained by a 
purely band picture, with intragap non-linear excitations suppressed by three-dimensional effects. 
However, recent experimental results on the clcctroabsorption (EA) and photoinduced absorption 
(PA) in substituted PPP by Lane et al. Q| are explained by the presence of non-linear excitations, 
such as singlet and triplet excitons, and charged polarons. 

The aim of this paper is to clarify the role and importance of the low-lying non-linear excitations 
in PPP by calculating its electronic structure and NLO properties in a realistic Hamiltonian. The EA 
spectrum compares favourably with recent experiments. We identify the key states which participate 
in the NLO processes. Moreover, by calculating the particle-hole separation of these states, we 
identify the band gap as the threshold state whose particle-hole separation increases linearly with 
oligomer size. This enables a rigorous determination of the band gap to be made. 

Recently, a two state molecular orbital model was introduced Q to describe the _Bi„ and Ag 
states of the phenyl based semiconductors. In the current paper we introduce a more thorough 
parameterisation of this model by fitting to improved exact Pariser-Parr-Pople model calculations of 
the molecular building blocks (i.e. benzene and biphenyl) 0|. This model is then solved for oligomers 
of arbitrary length without further parameterisation. 

As well as our earlier work, which was the first to use the DMRG method for the phenyl based 
semiconductors ^, there have been a number of other theoretical calculations on PPP. Bredas 
has used the VEH pseudopotential technique Champagne et al. have performed Hartree-Fock 
calculations and Ambrosch-Draxl et al. have performed density functional calculations using 

2 



LAPW and pseudopotentials [||. Rice et al. |jT^ have developed a phenomenological, microscopic 
model based on the molecular excitations of benzene. The absorption bands are calculated using an 
approximate Kubo formalism. In a series of papers, Shimoi and Abe have considered the optical and 
electro absorption of PPV in a re-parameterised P-P-P model using single-excitation configuration 
interaction (S-CI) ^ This method was used also by Harigaya et al. for studying optical 

absorption spectra and exciton properties in PPP, PPV, and related polymers |l3| . 

The structure of a PPP chain is shown schematically in Fig. 1. It possesses I?2/i symmetry. The 
electronic states can be classified according to their spatial, spin and particle-hole symmetries. In 
this paper we only consider states which are symmetric under refiection in the x — z plane, and 
either symmetric or anti-symmetric under refiection in the x — y plane. These will correspond to 
the low energy excitations. The ground state belongs to the (spin singlet, space- and particle- 
hole-symmetric) symmetry sector. Low-lying one-photon excitations occur between the and the 
^B^u (spin singlet, odd space and particle-hole symmetry) symmetry sectors. These excitations are 
polarised along the long (z) axis. We note that these states are electronically decoupled from the 
higher lying i?2u (symmetric under x — y refiection and anti-symmetric under x — z reflection) and 
i?3g (anti-symmetric under x — y and x — z reflection) states. 

Non-linear processes in polymers with inversion symmetry are determined by the third-order 
susceptibility x^^\~'^i — UJ2 — "-^a; "^i, ^2, ^3), which can be calculated as a sum over intermediate 
states (see, e.g., |ljl). It was suggested several years ago that most of the NLO properties can be 
described by an "essential states" model 17, The concept of the essential states model 

implies that there is a restricted set of states responsible for most of the NLO properties of the 
system. Usually, these are the ground state (1A+), the lowest excited (excitonic) state, the 

m^A+ state, which is the state most strongly related to the ^B^^ states through one-photon 
excitations, and the band threshold n^B^^ state. We investigate the validity of this approach for 
PPP. 

Until recently, numerical investigations of one-dimensional systems with strong electron-electron 
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interactions were limited to exact diagonalisations of short chains, to approximate configuration 
interaction calculations or to the use of density functional theory. A key aspect of this work is that 
we perform essentially exact calculations on our model using the density matrix renormalisation 
group (DMRG) method. We check the numerics by comparing DMRG results with exact results in 
the non-interacting limit, by comparing DMRG and exact diagonalisation calculations for a 7 unit 
oligomer and by monitoring the convergence of results with the parameter m (which controls the 
amount of Hilbert space truncation and hence the DMRG error) for longer systems in the interacting 
case. 

The plan of this paper is as follows. In the next section the molecular orbital model will be 
reviewed and parameterised. In §3 the energy and root-mean-square size of the low-lying states are 
calculated as a function of oligomer size, thus identifying the band threshold. The single particle 
excitation weight of the wavefunctions is also investigated. In §4 the dipole moments between states 
are calculated, thus identifying the key states in the NLO processes. The linear susceptibility is also 
calculated. In §5 we turn to the calculation of the EA and THG spectra using the sum-over-states 
method. Finally, we conclude in §6. 



2 The Molecular Orbital Model and Hamiltonian 

The starting point for the molecular orbital approach used in this paper is the well-known Pariser- 
Parr-Pople or extended Hubbard Hamiltonian: 

^ = - H *y [4cja + ll.C.] + C7^ - i j Ui^ -1)+IY1 ^'^("^ - - 

<ij>a i ^ ' ^ ^ ijij 

Here, cj^ and Cia are creation and destruction operators, respectively, for a tt electron with spin a 
on carbon site i, riic = cJ^Cjo- is the number operator, rij = nj| + riii and tij is the transfer integral. 
U and Vij are Coulomb repulsion parameters for electrons occupying one site and sites i and j, 
respectively. 
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The Ohno parameterisation for the Coulomb interaction is, 

= [i + ariyn' 

where a = (?7/14.397)^, thus ensuring that Vij e^/(47reorjj) as —^ oo, and is the inter- 
atomic distance in A. The optimal parameterisation for PPP, which was derived in 0, is [/ 10.06 
eV,the phenyl bond transfer integral, tp — 2.539 eV and the single bond transfer integral, tg = 2.22 
eV. 

2.1 The Model 

The six atomic 7r-orbitals of a phenyl ring cl^, i — 1, 6 may be transformed to six molecular orbitals 
(MOs), ajjg.. This approach for phenyl based semiconductors was used earlier by Soos et al. [|l9| . 
Recently, Chandross et al. po| employed a similar approach in their work on the characterization 
of excited states in conjugated polymers. It was suggested in |^ that the low lying Ag and Biu 
excitations can be described by only two of the resulting MOs, namely, the bonding eig HOMO 
and 6211 LUMO states. This is the so-called two state molecular orbital (2-MO) model. The other 
HOMO and LUMO states are non-bonding, because the wave function amplitude on the bridging 
carbon atoms is zero. Transitions between these states lead to high lying localised Biu excitations, 
while transitions which mix the bonding and non-bonding orbitals lead to excitations with B2u smd 
B^g symmetry We will not be concerned with these latter states in this paper. The (occupied) 
a2u and (empty) b2g states are situated far away in energy from the HOMO and LUMO, and are 
assumed to play only a minor role in the low-energy excitations. Below, the bonding HOMO orbital 
is denoted by |1) and the bonding LUMO orbital is denoted by |2). 

With two MOs left, and neglecting the three and four centre two-electron integrals and Coulomb 
interactions beyond nearest neighbour phenylene units, the transformation from the atomic to the 
molecular orbital basis results in the following Hamiltonian: 

H = - ^ ta/3 al^^ai+ipa + h.c.j + J2 ^ainia - 1) + C^Xl ("^"T - ^ j f'^^i " ^ j 

i a f3 a i a ia ^ 



l)("*+l/3 - 1) 





^ \ " t t 



(3) 



where S'^q = Epp' '^iap'^ pp'^iap' and <t are the Pauh spin matrices. 



The key interactions in this model are: the HOMO-LUMO gap (A = £2 — ei), direct onsite [U) 
and nearest neighbour {V) MO Coulomb repulsion, spin-exchange {X) and pair hop (P) between 
MOs on the same repeat unit, and hopping [t) between neighbouring repeat units. To understand the 
essential physics of this model, consider the limit t = 0. In this limit there are localised intra-phenyl 
particle-hole triplet and singlet excitations at %/ + p2 _ x and VA^ + + X, respectively. The 
MO Coulomb repulsion results in a potential well to the unbinding of these excitations. It costs an 
energy U ~V ~ X to seperate the particle-hole pair by one repeat unit and U — X to separate them 
by two or more repeat units. Finally, the hybridisation, i, leads to the delocalisation and ultimately 
unbinding of the particle-hole pair. 

A straightforward derivation of the new Hamiltonian parameters from the atomic Hamiltonian 
(1) gives results for the excitation energies which deviate from exact Pariser-Parr-Pople model cal- 
culations of benzene and biphenyl, as well as to over-estimating the optical gap by approximately 
1 eV in long oligophenylenes We therefore take the view that eqn. (3) contains the essential 
physics to model the low-lying excitations, but that these interactions are renormalised from their 
bare Pariser-Parr-Pople values. We parameterise the two state model by fitting its predictions to 
the exact Pariser-Parr-Pople model calculations of benzene and biphenyl . The two state model 
is exact in the limit that the interactions vanish, as in that case it describes particle-hole excitations 
from the valence to the conduction band. In the other extreme of the hybridisation vanishing, it 
correctly models localised intra-phenyl triplet and singlet excitations. In the intermediate regime the 
validity of the approach is determined by the test with experiment. We will show that the predictions 
are in good agreement with the experimental data, but in addition substantial physical insight is 
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achieved as to the nature of the excited states. 
2.2 Parameterising the Model 

Since the 2-MO model is apphcable to states of Ag and symmetry, it is parameterised by fits 
to states of the same symmetry in benzene and biphenyl. The interactions are parameterised in the 
following way: First, when the nearest neighbour hybridisation is switched off, the model should 
predict localised (phenylene) triplet and singlet excitons. The full Pariser-Parr-Pople calculation of 
benzene predicts a pair of excitations which are anti-symmetric under x — y reflection (the lEiu{z) 
and states) and a pair which are anti-symmetric under x — z reflection (the l_Bi„(y) and I-B2U 

states) in both the singlet and triplet channels Q. However, as was explained in ||2l|, the 2-MO 
model predicts that both the pair of lEiu{z) and 1-Bi„ states and the pair of li?i„(y) and I-B2M 
states are degenerate. In the 2-MO model the energies of the triplet and singlet are vA^KP^ — X 
and -y/A^ + p2 _|_ respectively. A and X (since P = X) are determined by fitting these values 
to the average values of the li?i„(z) and l-Bi^ triplet and singlet excitons obtained from the full 
Pariser-Parr-Pople calculation on benzene, which are 4.45 eV and 6.23 eV, respectively [Q. This gives 
A = 5.26 eV and X ~ 0.89 eV. Next, when the hybridisation is switched on the excitons dclocalise 
and interact. The key low-lying states are the long axis-polarised triplet [l^B^^) and singlet (l^Sj^) 
states, and the lowest even, covalent excitation, the 2^A+ [Q. Once A and X have been fixed, 
their energies are determined by t, U and V (we assume [tap] = t for all orbitals a, (3). We use the 
exact biphenyl calculations to fit these excitations. To simplify the fitting of these parameters we 
assume that V = U /2, and adjust t and U so that we have an exact fit to the full Pariser-Parr-Pople 
calculation of the biphenyl state and a minimum relative error for the l^B^^^ and 3^A^ states. 

This gives U = 3.67 eV, V = 1.835 eV, t = 0.895 eV, and a relative error for the 1^5+ and 3^^+ 
states of —2.2%. We note that since the band width is approximately equal to [/, these parameters 
are in the intermediate coupling regime. These parameters differ from those of Q, as they are based 
on a more rigorous fit to benzene and biphenyl. For ease of reference they are listed in table H. 
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2.3 Solving the Model 

(i) Density matrix renormalisation group solution 

The system is a one dimensional quantum lattice model with 16 states per repeat unit. For small 
lattice sizes, L, it is possible to calculate eigenvalues and eigenstates using exact diagonalisation. 
However, the largest system which can comfortably be reached by this method is the sexamer (L = 6 
phenylene units). In order to study larger systems, we turn to the DMRG method |2^. The 
DMRG is a powerful, robust, portable and highly accurate truncated basis scheme for the solution 
of low dimensional quantum lattice systems, and is especially well suited to the solution of open 
linear chains such as (j^). We have performed calculations of up to 10 low energy eigenvalues and 
eigenvectors in the ^A'^ and ^B^^ symmetry sectors, as well as dipole matrix elements, oscillator 
strengths and correlation functions for systems of up to 37 repeat units with sufficient accuracy to 
make comparisons with experiment. 

The DMRG is discussed at length in and reviewed in j24| so we restrict our discussion here 
to features relevant to our implementation of the method for (|^). The key features are the form 
of the system, environment and super blocks, the number m of states retained per block, and the 
good quantum numbers used to diagonalise the superblock Hamiltonian and the density matrix. 
We implement the DMRG for (|) using the infinite lattice algorithm |p3|| . That is, the system and 
environment blocks are reflections of one another, and arc increased by one repeat unit at a time, 
the initial blocks consisting of a single phenylene unit. The first two super blocks are schematically 
depicted in Fig. 2. They are comprised of the system and environment blocks abridged by a phenylene 
repeat unit. 

The total charge N — J2ia''^ia and the total z spin — ^ J2ia ('^mT ~ "•iai) are used as good 
quantum numbers in diagonalising the superblock hamiltonian and the system and environment 
block hamiltonians and density matrices. In addition, the spatial inversion (C2: aiaa a-L-i+iaa), 
particle-hole (J: aji^ 1-^ sgn{a)ai2g-, ajjo- sgn{a)aiia) and spin flip (P: Uiaa ^— > fliaff) symmetry 
operators can be constructed for the superblock. Their corresponding projection operators can 
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be applied to random superblock states in order to construct trial states of definite symmetries 
which can be fed into the sparse matrix diagonalisation routine used in the diagonalisation of the 
superblock hamiltonian. The resulting target states retain these symmetries as long as the iterated 
trial state is periodically resymmetrised, for example, every 30 or so matrix multiplications. This 
procedure is numerically stable because the density matrix eigenstates are eigenstates of the block 
symmetry operators and hence the obtained superblock energy eigenstates are exact (to within 
machine precision) eigenstates of superblock symmetry operators at each iteration. This procedure 
is checked for each calculated state by determining the expectation values of the symmetry operators, 
which are found to equal ±1 to around 12 decimal places. 

(ii) Accuracy tests 

We verify the validity of the DMRG solution by checking that the results obtained for the trimer 
and the pentamer agree with exact diagonalisation results. Basis truncation occurs for larger chains, 
the first being the septamer. This is the largest system that we can treat by exact diagonalisation 
(the dimension of the Hilbert space is 11 778 624). In Table ^ we compare the exact results with two 
DMRG calculations (with 65 000 and 130 000 states retained) for the energies and excitons sizes of 
the and 2^yl+ states of this system. It is clearly seen that, despite the fact that the DMRG 

calculation uses only a fraction of the total number of states (0.55% and 1.10%, respectively), in 
both cases the results are very close to the exact results, both for the energies and for the exciton 
sizes. 

In order to check the convergence for longer systems we first examinine the non-interacting 
(U — V — X = P = Q) case which can easily be diagonalised exactly for any chain length. In the 
DMRG calculations we retain m — 230 states per block. Exact and DMRG results are given in 
Table || for the ground and first excited state energies for a number of lattice sizes. We see that 
the DMRG resolves gaps between these states well and truly above the accuracy required in order 
to make comparisons with experiments, that is, a few hundreds of an eV. The accuracy is expected 
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to be even better in the interacting case where states are more locahsed and gaps are widened [£3j . 
That is, the accuracy should increase monotonicaUy as U is increased, up to the atomic hmit t^p = 
where exact results are trivially recovered by any real space renormalisation group procedure. In 
Table ^ we monitor the convergence of the optical gap in the interacting case with the truncation 
parameter m. We again see that the results have converged within errors which are negligible in 
so far as comparisons with experiment are concerned. Finally, DMRG calculations of the exciton 
correlation functions of the 1^-Bf„ and 2^A+ states are checked for systems of up to 21 repeat units 
by doubling the size of the supcrblock Hilbert space. The change in the exciton size is found to be 
less than 1% for l^B^^ and 3% for 2^A+ for all the systems studied. 

3 The Low Energy Spectra and Exciton Correlation Func- 
tions 

Table 4 shows the energies of the lowest triplet and singlet Biu excitons as well as the lowest covalent 
even C'A'^) excited state. The agreement between the predicted results from the 2-MO model and 
experiment are good for both oligomers and the polymeric limit, confirming the validity of our 
parameterisation. The calculated energies of the lowest excited and states as a function 

of oligomer size N are given in Fig. 3. Also, the charge gap defined as Eq = E{2N + 1) + E{2N — 
1) — 2E{2N), is plotted. Here, E{2N) is the ground state energy of a system with 2N electrons. The 
lowest excited state has ^B^^ symmetry, its energy always being lower than the charge gap. As the 
chain length is increased, the number of ^B^^ states below the charge gap also increases, creating 
an excitonic band in the limit iV ^ oo. (The states are inter-leaved with ^A^ states.) The 

energy of the lowest excited state, 2^j4+, almost coincides with the charge gap Eg, thus strongly 
implying that 2^A^ is a band threshold state. The energy of the lowest ^B^^ state above the 2^A+, 
denoted hereafter as n^B^^, also tends to Eq as iV ^ oo. A polynomial fit of the energies as a 
function of inverse oligomer length indicates that in the limit N ^ oo, the energies of the 2^A+ and 
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the n^B^^ states, and Eq all tend to a value of 4.04 eV. The energy of the l^B^^ state approaches 
3.30 eV. Thus, the energy results are good evidence that the lowest excited state is a threshold 
state dividing the spectrum into bound excitonic states below and unbound band-like states above 
it. 

That the S^A^ state is a band state is confirmed by an examination of the exciton spatial 
correlation function defined as: 

a,(|n)) = (n|4|li^+), (4) 

where 

Sjj = ;^(al2T«JiT + 42i"'3n) (5) 
is a singlet exciton creation operator, which removes a particle from the orbital |1) on site j and 
places it onto the orbital |2) on site i. Thus, 



gives the weight of single particle excitations in the state \ti) and 

p^j = cywi (7) 

is the distribution function for the particle- hole separation. Note that, as discussed in Appendix 
A, Cij = +Cji for states which are negative under the particle-hole transformation {i.e. 's'-wave 
excitons), while Ctj — —Cji for states which are positive under the particle-hole transformation [i.e. 
'p'-wave excitons). 

Using (7), we calculate the spatial extent of a given state, or the particle-hole separation, using 
the formula: 

l\\n))^{{^-Jf)^Y.^^-JfP,,. (8) 
(A similar approach to calculating the average particle- hole separation was used by D. Yaron and 



R. Silbey in their study of polyacetylene |28 .) The results are given in Fig. 4 as a function of 



oligomer size N . The ^B^^ states belonging to the excitonic band below the charge gap are those 
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with the smaUest electron-hole distance. With increasing oligomer size, the particle-hole distance in 
the l^B^^, as well as in other states in the excitonic band, approaches a limit of approximately one 
repeat unit. This indicates that these are strongly bound 's'-wave excitons. Conversely, the spatial 
extent of the 2^A+ and n}B^^ states are proportional to the oligomer size, and scales in the same 
way with N as do the lowest excited (unbound) states (namely, the l-Bj^ and 2A+ states) of the 
non-interacting model. 

We believe that all these results provide strong evidence that there exists an exciton in 

PPP, and that the lowest band state has the symmetry. The exciton binding energy for long 
oligomers approaches Eq — E^i^- — 0.74 eV. 

Finally, in Table 5 we show the single particle weight, Wi eqn (6), for the essential states (as 
defined in the next section) for a 15 site oligomer. Evidently, these states are predominately single 
particle in character. 



4 Oscillator Strengths and the Linear Susceptibility 

As a next step towards the calculation of the NLO properties of PPP, we compute the oscillator 
strengths of transitions between the essential states. These are, besides the ground state, the l^B^^, 
the 2^yl+, the n^B^^ and also the state situated above the n^B^^, which we denote as k^A^- 
This state has already been invoked for explaining the EA spectra of a number of luminescent and 



non-luminescent polymers by Liess et al. |29 . In Fig. 5 we show schematically the most important 
states for the non-linear optical properties. The results for the corresponding oscillator strengths are 
given in Table 6. 

Using the results for the oscillator strengths, the first-order optical absorption is computed. As 
expected, the dominant peak belongs to the lowest allowed transition, l^A^ l^B^^, while the 
band threshold state {n^B^^J is represented by a weak feature at 4.3 eV, as shown in Fig. 6 for a 
15 site oligomer. The higher lying optical transitions to states of i?2« symmetry are not indicated 
in this figure, as we do not model these states in our theory. The high lying localised 1^^^ — >^ B^^ 
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transition is also not explicitly indicated, but as discussed in §2, our theory predicts this to be at 



6.23 eV. 



5 Third Order Non-Linear Susceptibilities 

The NLO properties of PPP, such as third harmonic generation (THG) and electroabsorption (EA), 
can be related to the third-order macroscopic susceptibility x*^^'- The EA signal is related to the 
imaginary part of x^^'' (Hi- 

_Ar/r^l^W^)(-c.;c.,0,0), (9) 
cn 

while the THG spectrum is related to x*-^' (— 3aj; oj, uj). The third-order susceptibility x^^-* , in turn, 
results from the third-order microscopic hyperpolarizability ^xxxx'- 

X''xjxxi-^'y]t^l,^2,^^3) = -^fxx fxifxifxilxxxx, (10) 

where uJa- ~ +^^2 +^3 and the factor i results from the orientational averaging | |30t . The calculation 
of "ixxxx can be performed using the sum-over-states method (see, e.g., 



J / / ^J'gA^J■ABfJ'BC^^Cg ^ l^gAP-ABl^BCP-Cg 

^'^'^ AlBC ^^"^ ~ ^<y){^B - Wl - W2)(t^C - ^l) {^*A + ^i){.^B - UJi - L02){uJC " ^l) 
f^gAfJ-ABfJ'BCf^Cg fJ-gAf^ABfJ'BCf^Cg x 



-3 



(CJ^ + UJl){ujg +UJI+ UJ2){UJC - ^3) i^A + +UJl + ^2)(^^c + ^'^) 

^J■gA^J'Ag^J'gC^J'Cg . ^J■gA^^Ag ^^gC ^^Cg 



Yi- 

^ {U)A - CJ<t)(wa - t^3)i^C - t^l) (WA - W3)(cjJ + UJ2){UJC - t^l) 

^^gA^J'AgfJ■gCfJ■Cg , ^IgAtiAg^^gCKg f^-^^-^ 



+ uj^){ijj\+ uJ:i){Ljj}, + uji) {Ljj\+uj:i){ujc - ^2){(^c +^1)' 
where /i^j is the dipole matrix element for the transition between the states i and j, and K{—lJ(j; uji,lu2, W3) 
is a numerical constant which depends on the values of uJa, wi, uj2, 1^3 p^ - In particular, for the 
EA coefficient (ojo- = —lui = —lu,lu2 ~ UJ3 ~ 0), K = 3 and for the THG coefficient (ua- — — 3tJ,tJi = 
UJ2 = iL'3 — uj), K — ^. /i,2.3 denotes the average of all terms generated by permuting oji, UI2 and 
oj^. A finite lifetime of the levels A, B, C is taken into account in order to calculate ^xxxx at the 
resonance points properly. 
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The sum in equation (11) is over all states. However, due to the fact that the ground state belongs 
to the ^A^ symmetry sector, the dipole matrix elements are non-zero only for the transitions between 
and ^B^^. Thus, the states A and C in (11) are of ^B^^ symmetry, while the state B, as well 
as the ground state, are of ^A^ symmetry. 

Electroabsorption spectroscopy is used as a tool for investigating the low-lying electronic exci- 
tations in polymers (see, e.g., |3^, |3^, |l^). In |l^ it was shown that EA spectrum includes 
features of two types: First, those similar to the derivative of the first-order optical absorption; they 
are located at the energies of the levels. A second type are the peaks at the energies of the 
levels, arising from the spatial symmetry breaking in the applied electric field. Recent experimental 
studies of PPV and PPP EA spectra ^ ^ found both types of features in these systems. 

Using the transition matrix elements we calculate the EA spectrum of PPP oligomers using 
the 10 lowest eigenstates. The results for the 15 site oligomer are shown in Fig. 7. We note that 
transitions involving states of B2u and B^g symmetry are not included in this figure. The lowest- 
energy feature is a derivative-like feature arising from the ^-Bf„ level, showing a red shift. Next, the 
maximum corresponding to the 2^^+ state is situated at about 4.3 eV. The second derivative-like 
feature almost coinciding with the maximum can be related to n}B^^ state. Finally, the high energy 
feature at 4.8 eV corresponds to the k^A'^ state. The peaks corresponding to the 2^/1+ and ri^B^^ 
states are also identifiable in the EA spectra of Lane et al. at approximately 4.2 eV, just below 
the onset of the l^i?2« transitions. The high energy k^A'^ feature is not identifiable in their data, 
owing to the ^B2u transitions. For a comparison, results of the calculation with the five essential 
states only are also shown. Evidently, these states describe most of the EA spectrum. 

Finally, we show the third harmonic generation (THG) spectrum of PPP oligomers [Q. This 
spectrum is related to the x*''^H~3ci;, w, w, w). Mazumdar et al. used the THG spectrum in their 
discussion of the essential states [|l^. They argued that there are three main peaks in the THG 
spectra at E(liB{;)/3, ^{n^B^J/?, and ¥,{'m^A+)/2, where the m^A+ is the '^A+ state strongly 
coupled to the l^i?]T„ exciton. Our calculations (shown in Fig. 8 for a 15 site oligomer) confirm the 
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existence of these peaks (those at 1.15, 1.5, and 2.1 eV, respectively), the m^A+ state being the 2^A+ 
in our model. The peak at E(l^i?j^^)/3 is the highest in the spectrum. However, the calculation 
also reveals an additional peak corresponding to the k^A'^ state at 2.4 eV, as well as other peaks 
corresponding to ^A^ states. The feature at 3.4 eV corresponds to a one-photon transition to the 
l^Bj^ state. The essential states calculation shows features corresponding to the essentials states 
but fails to describe the total THG spectrum (employing 10 states) as well as it does for the EA 
spectrum. 

6 Conclusions 

This paper has been devoted to clarifying the nature of the low-lying A'^ and B^^ electronic excita- 
tions of PPP. A two state molecular-orbital (2-MO) model, whose parameters are derived by fitting 
its predictions to exact Pariser-Parr-Pople model calculations of benzene and biphenyl, was used. 
By solving this model using the DMRG method, it was shown that the lowest excited state is a ^B^^ 
('s'-wave) exciton, which is situated in the gap between the ground state and the band states. As 
the oligomer size increases, the number of excitonic levels also increases, thus creating an excitonic 
band in the limit N ^ oo. The calculated size of these excitons indicates that they are tightly 
bound, extending over only a few repeat units. The lowest band state is found to be a covalent 2^A+ 
state, whose energy almost coincides for all the oligomers studied with the charge gap Eq- Lying 
just above the 2^yl+ state is a band ^B^^ state (the n^B^^^ state). The size of the band states scales 
linearly with oligomer size. The binding energy of the l^-B^u exciton is thus determined rigorously as 
0.74 eV. The calculated values of the l^-Bj^ and 2^A+ energies of 3.3 and 4.0 eV, respectively, in the 
long chain limit, are in good agreement with the experimental results of 3.6 and 4.2 eV, respectively. 

The dipole matrix elements and oscillator strengths for the transitions between the lowest 
and ^-B^ states were calculated and the NLO properties of PPP, such as electroabsorption and third 
harmonic generation, were investigated. A comparison of the EA spectrum with the experimental 
data 1^ showed that the main features of the experimental spectrum are well described in the 2-MO 
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Hamiltonian. In particular, the derivative-like 1^-B{^ feature and the peaks corresponding to the 
band threshold 2^^+ and n^B^^ states are reproduced. 

Only five states account for most of the calculated electro-absorption. These are the , ^B^^, 
2^A+, n^B^^ and another band state, the A;^A+, thus confirming the essential states model. An 
analysis of the particle excitation weight of these states indicates that they are predominately single 
particle in character. 

In conclusion, the parameterised 2-MO model presented here gives a quantitative description 
of the low lying excitations and the NLO properties of the phenyl based semiconductors. The 
parameterisation was performed to achieve a good description of the states within the Ag and Biu 
symmetry sectors. Equally, however, by fitting to the relevant states of benzene and biphenyl, a 
parameterisation could be achieved for states in the B^g and B2u symmetry sectors. 
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A Particle-hole Symmetry in the Molecular Orbital Basis 

In this Appendix we derive the particle-hole inversion operator, J, for the molecular orbital (MO) 
representation, and discuss the particle-hole-symmetry-adapted wavefunctions for a chain. In the 
atomic orbital (AO) representation the action of the particle-hole inversion operator on the system 
can be cast as follows: 

Jc\^ = (-1)V,T, (12) 
Jc\^ = (-l)^+^ca. (13) 
Using the following relations between the MO and AO creation operators: 

-4^(2cL + cL-cL-24-cL + 4J, (14) 



12 

«L = ;^(2cL - 4. - cL + 24, - 4, - 4 J, (i5) 

we deduce the following formulae for the action of the particle-hole inversion on the MOs as: 

Ja\^ ^ a2i, Ja\^ = -a2'\, (16) 

= ai4, j4j^ = -ait- (17) 

Now let us consider the action of the particle-hole operator within a single exciton basis. Yaron 
and Silbey have considered the action of J on a single exciton basis using periodic boundary condi- 
tions Here, we develop these ideas using the real space MO representation. 

Let the ground state be represented by 

JV 

|G5)=n«liT4ilO). (18) 

1=1 

Then we may create a singlet excitation of spatial extent 5 — |i — i| and localised around the repeat 
unit n — {i + j)/2, as follows: 

€t7fl) = Sl\GS), (19) 



where 
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It is straightforward to demonstrate that 



J 



€-7/1) = -s], \GS) 



,n-<5/2\ 
€+S/2 ) 



(21) 



Thus, for an exciton locahsed about repeat unit n we may construct a state with definite particle-hole 
symmetry as, 



1^^) = E ( 



,n+5/2\ , 



' r 



.-S/2 
+<5/2 



(22) 



where the symmetric ('s'-wave) combination is negative under j, while the anti-symmetric ('p'-wave) 
combination is positive under J. It is = Cij which is essentially measured by eqn (5). 
A state of definite C2 symmetry may then be constructed as follows: 



W) 



/-J) sin ^ 



mrj 
N + 1 



(23) 
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Figure Captions 

Figure 1 : The structure of the poly(para-phenylene) chain. 

Figure 2 : Schematic representation of the superblocks used in the first two DMRG steps, showing the 
form of the system (S) and environment (E) blocks. 

Figure 3 : Energies of the lowest excited ^A^ and ^B^^ states relative to the ground state energy, and the 
charge gap Eg (dashed line) as functions of the oligomer size, N. lifii" (solid squares), 2^B^^ 
(solid diamonds), 3^B^ (solid triangles), 4^B^ (oblique crosses), n^B^^ (empty diamonds) 
and 2^A+ (empty squares). 

Figure 4 : Mean electron-hole distance for singly-excited low-energy states as a function of oligomer size, 
A'^. 1^-Bi„ (solid squares), 2^i?^ (solid diamonds), 3^5^^ (solid triangles), n^B^^ (empty 
diamonds) and 2^j4+ (empty squares). Results in the absence of electron-electron interactions 
are shown in dashed lines: l^B^^ (solid squares) and 2^^+ (empty squares). Also shown are 
dashed lines without symbols corresponding to N/2 and N/A. 

Figure 5 : Essential states for the non-linear properties of the PPP and the most important one-photon 
transitions between them. The oscillator strengths are shown for a 15 site oligomer. 

Figure 6 : Calculated first order optical absorption spectrum from the 1^^^ state to the low-lying ^B^^ 
states for a 15 site oligomer. 

Figure 7 : Calculated electroabsorption spectrum for a 15 site oligomer within the ^A^ and ^B^^ sym- 
metry sectors. Solid diamonds: full calculation, empty diamonds: essential states calculation 
(see text). 

Figure 8 : Calculated third harmonic generation spectrum for a 15 site oligomer within the ^^4+ and 
^Bj^ symmetry sectors. Solid diamonds: full calculation, empty diamonds: essential states 
calculation (see text). 
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Parameter 


Value (eV) 


A = 62 - ei 


5.26 


ill = ti2 = —t21 = —t22 


0.895 


X = P 


0.89 


U = 2V 


3.67 



Table 1: The values of the parameters used in the 2-MO model, eqn. (3). 



Property 


DMRG I 


DMRG II 


Exact 


Egs 


-26.821231 


-26.821232 


-26.821232 




-23.239216 


-23.239302 


-23.239318 


E2A1, 


-22.254639 


-22.254786 


-22.254843 


EiBi„ - Egs 


3.582015 


3.581930 


3.581914 


EiA^, — Eas 


4.560592 


4.5(i(i44(i 


4.5(i(i3<S9 


R.m.s.(lBi„) 


1.400599 


1.398925 


1.396940 


R.m.s.(2^g) 


2.660390 


2.652987 


2.653227 



Table 2: Energies of the ground and lowest excited states, as well as exciton sizes of the lowest 

excited states of a = 7 chain for exact and DMRG calculations (DMRG I: 65,000 states, DMRG 
II: 130,000 states). 
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N 


To (EXACT) 


To (DAiriG) 


Ey (EXiVC^T) 


Ei (DAiriG) 


3 


-16.9946597 


-16.9946597 


-13.9110457 


-13.9110457 


5 


-28.7708529 


-28.7708529 


-26.3070239 


-26.3070239 


7 


-40.5492152 


-40.5492152 


-38.3620419 


-38.3620023 


9 


-52.3277963 


-52.3277961 


-50.2911319 


-50.2910587 


11 


-64.1064055 


-64.1064050 


-62.1614635 


-62.1613673 


13 


-75.8850187 


-75.8850178 


-74.0003255 


-74.0002074 


15 


-87.6636326 


-87.6636311 


-85.8207064 


-85.8205619 


17 


-99.4422467 


-99.4422442 


-97.6294862 


-97.6293092 


19 


-111.220861 


-111.220857 


-109.430604 


-109.430389 



Table 3: Ground and first excited state energies in the non-interacting limit for various oligomer 
lengths A'' as calculated exactly and using the DMRG method. 



N 


m = 70 


m = 100 


m = 130 


m = 160 


m = 175 


m = 210 


7 


3.583 


3.582 


3.582 


3.581 


3.581 


3.581 


9 


3.480 


3.477 


3.477 


3.477 


3.477 


3.477 


11 


3.424 


3.420 


3.420 


3.420 


3.420 


3.420 


13 


3.392 


3.387 


3.386 


3.386 


3.386 


3.386 


15 


3.373 


3.365 


3.365 


3.365 


3.364 


3.364 


17 


3.361 


3.351 


3.351 


3.351 


3.350 


3.350 


19 


3.353 


3.341 


3.341 


3.341 


3.340 


3.340 
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?>:.m 


?,:.v.v^ 





Table 4: Values in eV of the exciton gap, the difference between the {l^B-^^) and (l^yl+) energies, 
as the truncation parameter m is increased, for a range of oligomer lengths N. 
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N 

1 V 


-Oil, 






H "V"T~* OT^I TT"1 On^" CI 1 1 IT^T 1 O 1 1 ^ CI 


9 


4:.OU 


fi 1 7 


o.oo 




o 




fin 


1 7 




A 


o.yo 


0. lo 


9 Q7 




o 


7Q 








R 

u 


fi7 


4 7n 


9 78 




7 


O.Oo 


4 p;7 


9 7*^ 




11 


3.42 


4.29 


2.65 




13 


3.39 


4.22 


2.63 




15 


3.36 


4.18 


2.62 




17 


3.35 


4.14 


2.61 




19 


3.34 


4.12 


2.60 




21 


3.33 


4.11 


2.60 




oo 


3.30 


4.04 


2.60 


3.43(b), 3.3(c), 3.5(d), 3.8(e) 



Table 5: Calculated vertical transition energies in eV for oligophenylenes of various lengths N . Note 
that the 2^A+ exciton is the lowest covalent Ag singlet excited state. Experimental results from 
biphenyl crystals (a) H and crystalhne films (b) |6|, (c) |3, (d) g, (e) §. 





nifii" (= 4iBrJ 


m^A+ (= 2^A+) 


k^A+ {= 4U+) 


0.996 


0.878 


0.794 


0.797 



Table 6: The single particle excitation weight, Wi (7), of the essential states for a 15 site oligomer. 
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N 


i'A+ ^ iisr„ 


liSfu ^ 2'A+ 


2^A+ ^ nifii" 


1^A+ ^ niSfu 


niSfu ^ k'A+ 


3 


21.330 (4.28) 


7.091 (1.32) 


26.391 (2.23) 


0.025 (7.83) 


21.202 (1.43) 


5 


35.858 (3.79) 


12.963 (1.10) 


48.603 (1.48) 


0.314 (6.37) 


30.152 (1.04) 


7 


49.237 (3.58) 


21.375 (0.99) 


66.291 (1.01) 


0.885 (5.58) 


40.231 (0.39) 


9 


62.077 (3.48) 


28.833 (0.92) 


79.324 (0.72) 


1.330 (5.12) 


77.925 (0.33) 


11 


74.603 (3.42) 


34.015 (0.87) 


87.975 (0.54) 


1.664 (4.83) 


107.969 (0.29) 


13 


87.109 (3.39) 


37.276 (0.84) 


93.555 (0.41) 


1.760 (4.64) 


124.494 (0.24) 


15 


99.662 (3.37) 


39.144 (0.82) 


97.084 (0.33) 


1.796 (4.51) 


133.758 (0.21) 


17 


112.411 (3.35) 


40.027 (0.80) 


99.091 (0.27) 


1.727 (4.42) 


140.365 (0.18) 


19 


125.153 (3.34) 


40.247 (0.79) 


100.176 (0.22) 


1.667 (4.36) 


144.788 (0.15) 


21 


138.013 (3.34) 


40.142 (0.78) 


100.612 (0.19) 


1.590 (4.31) 


146.118 (0.13) 



Table 7: Calculated oscillator strengths for selected transitions in oligophenylenes of various sizes, 



N. (The corresponding energy differences in eV are shown in brackets.) 
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